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ABSTRACT 

Wc consider warped equilibrium configurations for stellar and gaseous disks in the 
Keplerian force-field of a supermassive black hole, assuming that the self-gravity of the 
disk provides the only acting torques. Modeling the disk as a collection of concentric 
circular rings, and computing the torques in the non-linear regime, we show that 
stable, strongly warped precessing equilibria are possible. These solutions exist for a 
wide range of disk-to-black hole mass ratios Md/Mbh, can span large warp angles of 
up to ± ~ 120 deg, have inner and outer boundaries, and extend over a radial range 
of a factor of typically two to four. These equilibrium configurations obey a scaling 
relation such that in good approximation (f>/Q oc Md/Mbh where <f> is the (retrograde) 
precession frequency and is a characteristic orbital frequency in the disk. Stability 
was determined using linear perturbation theory and, in a few cases, confirmed by 
numerical integration of the equations of motion. Most of the precessing equilibria 
are found to be stable, but some are unstable. The main result of this study is that 
highly warped disks near black holes can persist for long times without any persistent 
forcing other than by their self-gravity. The possible relevance of this to galactic nuclei 
is briefly discussed. 
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1 INTRODUCTION 

The increasing power and spatial resolution of modern 
observations has provided evidence that warps are not 
unique to galactic disks, but appear also on much smaller 
scales. These include nuclear and accretion disks surround- 
ing supermassive black holes in galactic nuclei (nuclear 
disks hereafter). The pioneering example is the maser 
disk of NGC4258. T he high velocity masers mapped by 
iMivoshi et all (|l995h are best explained by the existence 
of a mildly warped disk extending from 0.13 to 0.26 
pc |Herrnstein et all Il996t ). The nearby Seyfert galaxies 
NGC10 68 and Circinus also harbor warped disks i n their 
centers (|Greenhill et al ] |2003al : iGalhmore etaHl2004 ) , again 
traced by the maser emission. Of the ~ 100 massive young 
stars in the center of our Galaxy about half form a well- 
defined, warped disk, and some of the others are on a 
counterrotating structure which may be a dissolvin g disk 
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Nuclear disks can develop a warped shape through sev- 
eral mechanisms. Very close to the cente r, the dragging of 
inertial frames by a rotating black hole (|Lense fc Thirrind 
1 1918h causes precession of a planar disk, if it is inclined to 
the plane perpendicular to the black hole's spin. Internal vis- 
cous torques try to align the disk angular momentum with 
the black hole angular momentum. Beyond a transition ra- 
dius, the disk does not feel the effect of the black hole and 
remains at its initial inclination, while inside this radius the 
alignment proceeds. Hence the dis k becomes warped [the 
iBardeen fc Pettersonl (|l975l ) effect] . iNataraian fc Armitagd 
jl999l ) showed that for black holes with masses of order 
10 s Mq and accretion rates close to the Eddington limit 
the alignment time scale is short (t ^ 10 6 yr). Application of 
this effect to the warped disks of NGC4258 and NGC1068 
shows that the alignment radius lies well inside the observed 
positions of the maser spots, and model s can be constructed 
that fit the o bserved warps rather well |Caproni et alj|2007l : 
lMartinll200d h 

When a warped disk is exposed to radiation from a cen- 
tral source, or from its own inner portions, it is not illumi- 
nated isotropically. If it is also optically thick, the emission 
received at each position and re-radiated perpendicular to 
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the local disk pla ne induces a tor que on the disk, and the 
warp is modified (|Pettersonlll977l ). Perturbations t o planar 
disks can therefore cause radiation driven warping (jPringld 
1 19961 ). Assuming a radiative efficiency e ~ 10 2 , and a black 
hole mass of 10 8 Mq, an initially flat disk is prone to warp- 
ing beyond r ^ 0.1 pc, when the vertical an d radial viscosity 
coeffic ients are comparable (|Pringlei[l997h . iMalonev et al.l 
|l996l ) studied the stable and unstable modes of radiatively 
excited linear warps and found that the warp in NGC4258 
may be explained by this mechanism only if the radiative 
efficieny is high. 

Warps generated by gravitational interactions 
have been invest i gated mainly in the galactic context. 
iHunter fc Toomrel |l969) studied the linear bending waves 
of a self-gravitating, isolated, thin disk. They showed 
that such a disk permits long-lasting bending modes only 
when its surface density near the outer radius is truncated 
sufficiently fast, but not when realistic smooth edges 
are considered. This suggested interactions with nearby 
companion galaxies as a likely cause of warp excitation. 
Later, as the evidence for dark matter halos around galaxies 
became stronger, modelers developed scenarios in which the 
disk assumes the sha pe of a norma l mode in the p otential of 
a flattened dark halo (Sparkc 1984; Kuiikcn 1991). However, 
subsequent work showed that these modes are damped 
quickly when the in ternal dynamical response of the halo is 
taken into account jNelson fc Tremain el ll995l ; iRrnnev et al.l 
1 19981 ). Today, it seems most plausible that galactic warps 
result from interactions and fro m accretion of materia l with 
misaligned angular momentum (| Jiang fc Binnevlll999l ). 



On nuclear disk scales, iPapaloizou et afl (|l998i) stud- 
ied in linear theory the evolution of a thin self-gravitating 
viscous disk interacting with a massive object orbiting the 
central mass, with application to NGC4258. They concluded 
that the warp in the maser disk of NGC4258 might have 
been excited by a binary companion with a mass compara- 
ble to or higher than that of the maser disk. Their model 
also sugg ests a small twist (i.e. varying line-of-nodes) due to 
viscosity. iNavakshinl l|2005h considered the case of a non-self- 
gravitating disk perturbed by a massive ring. Employing the 
gravitational torques in the linear regime, he evaluated the 
precession induced by the ring on the disk elements. When 
the self-gravity of the disk is not taken into account, the 
rings precess differentially, which tends to destroy the disk 
structure. 

Can models of warped nuclear disks be generalized to 
the fully non-linear regime? And assuming that the observed 
warps in galactic nuclei have been excited by one of the 
mechanisms discussed above, can the disk self-gravity main- 
tain the warp even after the exciting torque has ceased to 
exist? As a first step towards answering these questions, 
the goal of the present paper is to investigate the possibil- 
ity of steadily precessing, stable, non-linearly warped self- 
gravitating disks in the (Keplerian) gravitational potential 
of a massive black hole. In the following sections, we use a 
simple circular orbit ring model to find stable warped equi- 
libria for systems with 2, 3, and many rings, assuming that 
the self-gravity of the rings provides the only acting torques. 



2 STEADILY PRECESSING WARPED DISKS 
AND THEIR SCALING RELATION 

2.1 Cold Disk Model and Equations of Motion 

We consider a cold disk in which stars or gas are assumed to 
move on very nearly circular orbits. Following similar anal- 
ysis of galactic warps (e.g., Toomre 1983, Sparke 1984, Kui- 
jken 1991) we model such a disk as a collection of concentric 
circular rings. The orbital motion in the disk is maintained 
by the central black hole, and the self-gravity of the disk 
causes the rings to precess around the total angular mo- 
mentum direction. Each ring may represent a set of stars 
or gas elements uniformly spread around their circular or- 
bit. Moreover, when the precession frequency arising from 
the self-gravity of the disk is small compared to the orbital 
frequency of motion, the orbital parameters of single stars 
change only slowly and so one can average over the orbital 
motion. In this case, also the force exerted by a single star 
or mass element on the rest of the disk can be replaced by 
the force due to a ring of material spread over the orbit 
(|Goldreichlll966h . 

Any of the rings is characterized by its mass mi, radius 
n, inclination angle 9i with respect to the reference plane, 
and azimuthal angle <j>i where the line-of-nodes cuts this 
plane. Later we will identify the reference plane as the plane 
perpendicular to the total angular momentum vector. T he 
Lagrangian d of ring i is given by (|Goldstein et alj|2002h : 



d = ^ 2 +, 



>i sin 2 0i) + m ' r ' (i> + <}>i cos 9i) 2 



V(r 



,4>i 



(1) 



The first two terms in equation ([T]) represent the kinetic 
energy of the motion Ti, V(ri, 0i, 4>i) represents the gravi- 
tational potential energy, the Lagrangian is d — T — Vi, 
and the energy of a ring is Ei — T + Vi . tp is the position 
of a point on the ring, measured from the ascending node; 
(9, <p, ip) are Euler angles. The angular momentum of the 
motion along the ring 



= m 4 r i n(ri) = m 4 r 2 (i/; + 4>iCos8i) 



(2) 



is conserved since d does not depend on the coordinate 
tp. The other momenta are the p ( f, i , the angular momentum 
around the z— direction, and pg i , the angular momentum 
around the line of nodes. The equations of motion are: 



P<t>i 



rrurf 



Hi i : . 2 a i , 
— — (pi sin 0i + p^ i cos I 



(3) 
(4) 



rmri ; 2 



sin t>i cos t>i — q>iP^, i sin ( 



dV{ri,6i,<t>i) 



d9 t 



(5a) 



dV(n,9i,(pi) (p t j >i - p^, i cos9i)(p,p i - p<t >i cosOi) 



dV{nAAi 



(5b) 
(6) 
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and the Hamiltonian is: 



Hi 



2.2 



P% 1 



P%i 



rrurf 2 rmr'. 



+ 



(jPfc -Vi H cos 8 iY 
rmr 2 sin 2 Oi 



+ V{r i ,e i ,4> i ). 



(7) 



Components of V(r, 8, < 
the Torques 



and Evaluation of 



The gravitational potential energy, V(r, 8, <p), has two com- 
ponents. One arises due to the central black hole, and is 
simply 

GmiM bh 



V bh 



r, 



(8) 



at the position of the ring. The other component is the po- 
tential term V m arising from the interaction of the ring under 
conside ration with all other rings . We follow the descrip- 
tion of | Arnaboldi fc Sparkd |l994j), using the derivation of 
iBinnev fc Tremainel (|l987l ') (Section 2.6.2), to evaluate the 
torque arising from the ring interactions. 

The gravitational potential due to a circular ring of 
mass mi and radius in the (x, y) plane is 



$>(x,y,z) 



2Gm, K(k)^(l-k 2 /2) 



where 



fc 2 = 



4Rr, 



(r 2 +r 2 +2R ri ) 



(9) 



(10) 



Here K(k) is the complete elliptic integral of the first kind, 
and R is the cylindrical radius R 2 — x 2 + y 2 , so that R 2 = 
r 2 — z 2 . A second ring of radius r, at an angle onj to the 
first ring follows a curve z = r 3 sin aij sin ip, where tp runs 
between and 2tt. The mutual potential energy is 

V^) = - 7r2( ^^ )1/2 fj K{k)^W]2d^ (11) 

where raj is the mass of the second ring, and fc depends on 
ri/rj, sin and ip. The angle aij between the two rings is 
given by 



cos aij = cos 8i cos 8j 

which reduces to cos aij 
nodes are aligned {(pi = 
rings is 



f sin 8i sin 8j ■ cos(<^ — <j>j), (12) 

= cos(#i — 8j) when the line-of- 
fij). The torque between the two 



dVij 
dai. 



GmirajTiTj sin 2atj 
(r\ + r 2 )3/ 2 



Iij {aij , i~i/rj), 



(13a) 



ir/2 



E{k)(l - k 2 /2) 



(1 



(1 

fc 2 /2) 3 / 2 



fc 2 ) 



K(k) 



sin 2 tp dip 



k 2 n ~ r^y- (13b) 

K y 1 — sin atj sin ip 

We u se the numerical program of lArnaboldi fc Sparkd 
(| 19941 ) for evaluating the integrals in this expression. The 
torques with respect to the angles (8i, (pP) follow from multi- 
plying equation (|13a|l by daij /ddi or daij/dcpi. In the follow- 
ing, we will write V m ,i = Ylj=ti l° r the potential energy of 
ring i due to the other rings, so that its totaJ potential energy 
becomes V{ri,8i,cpi) = Vbh{ri) + V m ,i- For further reference 
we also define My = —dVij/ddi, and Mg,% = —dV m ,i/d8i 
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Figure 1. The mutual torque between two neighboring concentric 
rings, whose radii r ou t and r; n are in the ratio v as specified on the 
plot. The linear regime is limited to the left of the peak amplitude 
in these curves. 



for the total gravitational torque on ring i around its line- 
of-nodes. 

Figure [T] shows the torque between two rings with radii 
in the ratio v = r out /ri n as a function of their mutual incli- 
nation a. The maximum of the torque occurs at very small 
angles, as noted previously by Kuijken (1991) who gives the 
approximation a max — 1.2\u — 1|. Only for a < a m ax can 
the mutual torque be approximated as a linear function of 
a. Thus solutions 8(r) for the warp shape in linear theory 
can be scaled by a constant multiplicative factor only so long 
as the local gradient dd/dr < 1.2/ r. Otherwise the local self- 
gravity torques of the disk are no longer able to maintain 
the linear theory warp shape, the non-linear equations must 
be used, and the shape of the warp must change. 



2.3 Steadily Precessing Equilibria 

A configuration of inclined rings precessing as a rigid body 
with constant (p will in the following be denoted as a steadily 
precessing e quilibrium, or equi l ibrium for sh o rt. In ear- 
lier works bvlrlunter fc Toomrd (|l969T l; ISparkd (|l984h . and 
ISparke fc Casertand (|l988l ) it was found that the eigen- 
frequencies of linear m =l warp modes are purely real. 
IPapaloizou et all (|l998l ) showed that this is a consequence of 
the self-adjointness of the operator in the tilt equation when 
there are no viscous or non-conservative forces. Then the 
eigenvectors are also real and thus the warp has a straight 
line-of-nodes and no spiralitjQ. In the light of these linear 
theory results our effort will also be to find equilibria where 
all the rings have the same azimuth (p. The condition that 
all rings maintain constant inclination, 8 — 0, implies also 
pe = 0, and for a given precession rate (p, the simultaneous 
solution of this equation for each of the rings determines the 
inclination angles, i.e. the equilibrium shape corresponding 



The line-of-nodes for a set of orbits is here defined as the 
union of points where all inclined orbits (0; ^ 0) with respect to 
the reference plane intersect this plane. 
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to this value of <j>. We note that for <fi = equations (|5a|l ad- 
mit a trivial tilt solution 8i — const, but will assume <j> 7^ 
in what follows. 

From equations (|5a,|) and Q we can solve for the pre- 
cession rate of ring i: 



COS 8; 



+ 



rriirt cos ( 



(14) 



when 9i ^ 0. Here f2; = y GM^h / r\ is the angular velocity 
of particles on the ring around the black hole, and the term 
Ylj dVij I d9i is the torque on ring i caused by all other rings 
j. The precession rate can therefore be fast or slow, cor- 
responding to the plus and minus signs in this expression. 
When the interaction potential Vij increases away from the 
plane 8 — 0, the second term in the square root is positive, 
so that the slow precession is retrograde (<f) < 0). In the 
remainder of this paper we focus on such slow retrograde 
precession. 

The components of angular momentum along the orig- 
inal (x, y, z) axes for a single ring read: 



Ixi = PBi cos <f>i + 



sin 0i(p^ i - V<t>i 



pe i sin < 



cos cf>i (p^, i — p^ cos 8i 



sin Oi 



1^ = P<t>i 



(15) 



(16) 



(17) 



Let us assume that we have found a precessing equilibrium 
from solving equations I] 316(1 . with pe i = 0, pe i = 0, and 
p ( p i = const., (j> = const., cj>i — (f>. Inserting equations (|5a[) 
and ((4]) into the expression for l Xi , simplifying and summing 
over all rings gives the total angular momentum 



■ V —v 







(18) 



which sums to zero because for each pair of rings with in- 
teraction potential Vij the torques are equal and opposite. 
Similarly, the total l y = 0. Thus the total angular momen- 
tum of such a precessing equilibrium configuration is paral- 
lel to the z-axis. By construction, the angular momentum of 
the precession alone is also along the z-axis, i.e., the disk 
precesses around the total angular momentum vector axis. 

For a uniformly precessing configuration, additional in- 
sight may be obtained by moving to a coordinate system 
which rotates around the angular momentum axis with the 
disk's precession frequency cj> (Kuijken 1991). In this refer- 
ence frame the shape of the precessing disk is stationary, 
but the particles in the different rings still spin about their 
rings' symmetry axes. If the particles in ring i rotate with 
velocity Q.(ri)ri in the positive sense, they experience a Cori- 
olis force in the rotating system which, integrated over the 
ring, results in a Coriolis torque on ring i along the pg-axis 
(line-of-nodes), given by 



Mc,i = — m i r i n(r i )(/) sin ( 



(19) 



For < 8 < 7r/2 and negative <j> this torque is along the 
positive pe-axis, i.e., is trying to retard the ring relative to 
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Figure 2. Relative inclination of two rings precessing together 
around a central mass M^, as a function of the square root of 
the ratio of their radii. The combined mass of the two rings is 0.5% 
of the central mass, approximately as in the Galactic center. The 
different curves are for various values of the precession frequency 
of the two-ring system, in units of the Keplerian frequency at r\ . 



the rotating frame. Because the retrograde precession speeds 
are small, we can neglect the centrifugal force terms. In this 
case, a stationary precessing configuration is obtained when 
the forward gravity torques and the retarding Coriolis torque 
balance in the rotating frame. 



2.4 2-Ring and 3-Ring Cases 

The argument just described suggests that there should exist 
steadily precessing 2-ring configurations in which one ring 
is tilted above the plane 8 — and a second ring is tilted 
below this plane. Both rings are pulled towards 8 — by 
the gravitational force from the other ring. The resulting 
gravity torques cause the angular momentum vectors of the 
two rings to precess in the same sense, and are balanced by 
the Coriolis torques in the precessing frame. To find such 
configurations we need to solve pg = using eq. (|5a|l for 
both rings simultaneously. Assuming <j> <C O, we can neglect 
terms of order 2 ; then using eq. ((2} the equation for the 
inner ring at r\ becomes 

sin 6>i 



-QQ-/mir 1 il(n) 



and the ratio of the two equations is 
sin 8 1 j sin 82 



m2rlQ{r2)/m 1 rlQ,{r 1 ) 



(20) 



(21) 



where mi, 7712 are the two ring masses, n, T2 their radii, 8\, 
82 their inclinations, and V12 the interaction potential. Given 
the ring masses and radii and 81, say, we can determine from 
these equations 82, the interaction potential, and thus finally 
the precession rate <j) required for steady precession. 

Using the expression in (|13a|l for the torque between 
the rings, equation (|20|l can be cast into a more useful form: 



sin 2a 7i2(a, v) 



sin 0i (1 + ;/ 2 ) 3 / 2 M bh ' [ ' 

where the angle a = 9x — 62, v = rz/ri, fi = ma/mi, and 
7i2 denotes the integral expression of equation (I13bl l. 

Figure [2] shows the difference 8\ — 82 between the incli- 
nation angles of the two rings versus the square root of the 
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Figure 3. 3D view of a 3-ring system. The middle ring lies close to 
the equator, while the others are distributed almost symmetrically 
around it. 



ratio of their radii, Ti/r?, for different precession frequen- 
cies, expressed in units of the Keplerian frequency at n. 
The combined mass of the two rings is chosen to be 0.5% of 
the central mass, approximately as infer red for the system o f 
two stellar rings in the Galactic center (|Genzel et al.|[20"03h . 
6\ — 62 increases with decreasing precession speed when the 
mass ratio is fixed. Eq. (|22|) shows that the same process- 
ing equilibrium configuration can be obtained by changing 
(j> oc Ma = nil + ni2 and leaving all other parameters un- 
changed. More massive rings must precess faster for the same 
inclinations. Thus the sequence of curves in Fig. [2] can also 
be interpreted as a sequence of fixed precession frequency 
but with mass ratio Md/Mbh increasing from bottom right 
to upper left. 

Next consider three rings. In this case, each of the rings 
precesses in the potential of the other two rings, and the 
reference frame is defined by the common plane of precession 
perpendicular to the total angular momentum vector. Again 
pe [eq. (|5a|l ] should be zero at equilibrium for each of the 
rings. We can sum these three equations: 



3 3 / 2 

^Pe s = I ' 1 4? sin 8i cos di - <j>p^ t sin di 



1=1 jy^i 



(23) 



The Vij terms cancel since dVij/dOi = —dVij/ddj. The re- 
maining terms can be rewritten as 



y rriirl sin di 



-tl(ri) + -0cos#i 



0, (24) 



making use of eq. This shows that, apart from the no- 
precession solution, a steadily precessing equilibrium is pos- 
sible only when at least one of the rings lies on the op- 
posite side of the equator with respect to the others, i.e., 
has 9i < 0. Likewise the two rings of a precessing two-ring 
system must lie on opposite sides of the equator. Figure [3] 
shows as an example the 3D view of a 3-ring system with 
mass M d = 0.05Af Wl and <j> = -0.0021ft(r 2 ). 
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Figure 4. Inclination of a disk of constant surface density at 
different radii. The model consists of 35 rings processing at a rate 
of rfi = — 10~ 3 units, so that <p/fl(r Te {) = —0.0021 on its middle 
(reference) ring. Each curve corresponds to a different Md/M^u 
mass fraction. The warp becomes more pronounced when the disk 
mass is increased. The smallest and highest masses correspond to 
the limits for stability (see Section 12.71 1. 



2.5 Approximation of A Disk With n-Rings 

We now consider a disk represented as a collection of n 
concentric rings. To find a precessing equilibrium, we solve 
pe = (eqs. I5a[) for all rings simultaneously, summing over 
the torques from all other rings (eas. ll3a")l . These are n equa- 
tions for n+1 unknowns, the n inclinations dj and (f> . which 
we solve keeping cj> fixed (|Arnaboldi fc Sparke|[l9 94 rl. Fig- 
ure [4] shows a sequence of equilibria obtained for a constant 
surface density disk consisting of 35 rings. On each curve, 
the extent of the disk (i.e. Ar — r out — ri n ) is fixed at 2.2 
units, and the precession rate is 0/f2(r ro f ) = —0.0021 where 
f2(frof) is the circular frequency on the middle (reference) 
ring. The disk mass fraction Md/Mbh varies from 0.16% to 
21%. As the mass of the disk increases, the degree of warp- 
ing increases dramatically so that the Coriolis torques can 
keep the balance of the gravity torques. The basic shape of 
the disk is similar to that of the system of three rings in 
Fig. [3] The middle rings lie closest to the equator, while the 
inner and outer rings are almost symmetrically distributed 
around it. 

Obviously, the larger the number of rings the better 
the approximation to a continuous disk. Figure [S] shows the 
convergence of the total torques (upper panel), and of the in- 
clination angles obtained in steadily precessing equilibrium 
(lower panel), for the innermost and outermost rings, when 
the number of rings to represent the disk is increased but 
the extent and the mass of the disk are kept fixed. One sees 
that quite a number of rings are needed before the torques 
converge. The inclination of the outer and inner rings have 
approximately converged when n > 30. 



2 Note that this does not work in linear theory because the lin- 
ear solution can be scaled arbitrarily, i.e., one of the di can be 
eliminated. 
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Figure 5. Convergence of the total torques (upper panel), and 
inclinations (lower panel), for the innermost (solid lines) and out- 
ermost rings (dot-dashed lines) , when the number of rings to rep- 
resent the disk is increased but the extent and the mass of the 
disk are kept fixed. 



2.6 Scaling the Solutions 

Now we go back to the equilibria themselves, in particular 
to the question of their scaling properties. When the torque 
on ring i from all other rings is decomposed as 

dVij 



M G ,i 



EdVij dctij 



(25) 



• precession rate - m nr 




the mass- and radius-dependent part is the first derivative 
on the r.h.s. The second factor in each term of this sum 
depends from equation (|12|l only on the two sets of angles 
9i , <j>i , 9j , <j)j . For equilibria with a common precessing line- 
of-nodes, ay — (9i — 9j), so the derivative is always unity. 
For the potential derivative terms, we use equation (|13a[) 
and express all ring masses and radii in terms of the mass 
and radius of a reference ring, i.e., we write [ii = mi/m Te {, 
Vi = ri/r Ic f, and similar for j. Then (|25[) takes the form: 



M G ,i = 



Gm 2 ci 
r lc f 



E 



HiUji/ii/j sin oiij 



(26) 



where we denote the expression over the brace as Di/2. For 
each precessing equilibrium disk configuration as in Fig. [4] 
Di is a constant, and the torques on all rings scale as tx 
rn ref/ r ref ■ If we now go back to equation (|14[) . insert equation 
()26[) . and normalize the precession rate with the circular 
frequency at the reference radius, Q rc f = (GMbh/r^f) 1 ^ 2 , 
we find 



3/2 
V- COS Oi 



1± ,/l 



Vi COS Oi m Te f 

~Mbh 



fj,i sm ( 



(27) 



For negative (j> and after a Taylor expansion of the term in 
the square root, appropriate for slow retrograde precession, 
equation (|27|) becomes 



Di 



fir 



m rcf 



2^!A L/2 sinf3 i M bh 



Di m Tci 
2Ai M bh 



(28) 



A precessing equilibrium is one for which all rings precess 
with the same common frequency, <f>i = cj>- Equation (|28[) 



Figure 6. Top: Values —Di/2Ai (eg. I28| l for a system of 5 rings 
with the parameters described in the text. The absolute value 
of the precession speed increases along the curves from top to 
bottom. Bottom: Parameter scaling of this system according to 
cq. ||28|I . This figure shows a zoom into the region where the 
orientation of the 1 st ring predicted by the scaling deviates from 
the orientation in the original 5 ring model (shown as the solid 
line) .Dots show when the black hole mass is increased by a factor 
of 4 and the ring mass is increased by a factor of 2, the dottcd- 
dashed line shows when the radius is increased by a factor of 4 
and ring mass by a factor of 8, and the long dashed line shows 
when the scaling is done by increasing the precession rate by a 
factor of 2 while increasing again the ring mass by a factor of 2. 



thus shows that for a fixed precessing disk mass configura- 
tion (i.e., fixed ring masses, radii, and inclinations, hence 
fixed Di/2Ai), the precession rate 4> scales proportional to 
the Keplerian frequency at some reference radius in the disk 
and proportional to the disk-to-black hole mass ratio. Vice 
versa, equation (|28|l can be interpreted as a scaling relation 
which says that a precessing equilibrium solution remains 
unchanged in shape ((?;) under changes of the disk mass, 
disk radius, black hole mass, and precession rate, provided 
the ratio (<j)/il Te t)/{Md/Mbh) is held constant. 

Figure [6] depicts the values of —Di/2Ai for a system of 
5 rings. The radii of the rings are calculated such that r; = 
x n, with i = 1, 2, ..n, K = 1.07, n = 5.75, and n — 5, 



so if the third ring is the reference ring, i>i = 



The ring 



masses are assumed to all have the same value, 0.5516, so 
Hi = 1, the black hole has a mass of 51.16, and Md/M^n = 
0.054. On each curve in Fig. [6] the precession rate <p increases 
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with steps of-5xlCT 5 (A0/fi(r 3 ) = -1.18xlO~ 4 ) starting 
from a value of -1 x 10~ 4 (i.e., 0/fi(r 3 ) = -2.4 x 10~ 4 ) 
at the top. We checked the accuracy of the scaling and of 
our calculations by computing —Di/2Ai values for different 
parameter pairs of the system that should give the same 
—Di/2Ai according to equation ()28p . We overlay the results 
for the first ring and precession speed 0/O(r 3 ) = —0.0021, 
in the lower panel of Fig. [6] zooming into the parameter 
region —2.3 < D\/2A\ < —1.9 where the different curves 
deviate from each other the most. In the worst case, due to 
the scaling of the ring radii, the deviations of the 8i's from 
their values for the original 5 ring system are still less than 
1°. Changes in radii cause the largest deviations from the 
scaling relation because of the way in which they enter in 
the quantity Di (equation I27p . The scaling results for the 
other rings are similar. 



2.7 Stability 

In this section we investigate the stabil ity of the process- 
ing eq uilibrium solutions found above. Iflunter fc Toomra 
(l969|) proved that isolated thin self-gravitating disks are 
stable to all m = 1 warp perturbations and this carries 
over to disks embedded in sp herical or oblate potentials 
fe.g. lSparke fc Casertandll988h . We show here that the non- 
linearly warped precessing disks can be both stable and un- 
stable to general ring-like perturbations. We describe small 
perturbations of the precessing disk solutions by the lin- 
earized equations of motion around equilibrium: 

d 2 T> 



Mi 
Apg t 
Ap^ 



dp 2 e . 
d 2 Tj 

d 2 T, 



d 2 T 
dp 2 



d 2 V m 

d 2 V m ,i 
d<j)id6i 



-Ap^ 



V4h 
d 2 T, 



A9i - 



d0 2 



Afa - 



d 2 V„ 



AO, - 



j 

d 2 V mA 



94> 2 



A(j)i 



Ed 2 Vm,i 



A0, 



E9 2 V m ,i 



Ah 



(29) 



Here Tj and V m ,i are the kinetic and potential energy terms 
in the Hamiltonian (O for ring i, respectively, and the par- 
tial derivatives are evaluated at the equilibrium solution 
(di, (pi = const., =0,p t f, i ~0). These linear equations have 
solutions of the form e At A#o, • ■ • etc., where A = \-n ± \\i m 
with its real and imaginary parts. The (constant) coefficients 
of the A terms in equation (|29|l form a matrix H which car- 
ries the information on stability. When the real parts of the 
eigenvalues of the matrix H, Xn = 0, the imaginary parts of 
the eigenvalues, \i m , constitute a rotation matrix through 
which the solutions oscillate around the precessing equilib- 
rium with frequencies Aim, and the equilibrium is said to be 
stable. When X-jz < 0, the solutions spiral towards the unper- 
turbed equilibrium positions, leading to asymptotic stability 
of the equilibrium. If, however, any of the eigenvalues have 
a nonzero real part, \n > 0, the system moves away from 
equilibrium exponentially, and is unstable. 
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Figure 7. Change of the shape of a solution near the lower stable 
boundary of Fig. [4] with the number of rings used to represent 
the disk. The disk-to-black hole mass ratio, precession speed, and 
number of rings are given on the plot. The letters "u" and "s" in 
parentheses following the number of rings stand for unstable and 
stable solutions respectively. 



For determining the stability of any of our precessing 
n-ring solutions, we compute the An x 4n stability matrix 
H, using the equilibrium (#;, </>;= const., pg i = Q,p ( j, i = Q). We 
then evaluate the eigenvalues of the matrix H, using routine 
F02EBF of the Numerical Algorithms Group (NAG). This 
routine is suitable for computing eigenvalues and optionally 
eigenvectors of real matrices. 

First, we briefly discuss one example for the conver- 
gence of the linear stability results. In Section[23]we had al- 
ready discussed the convergence of the gravitational torques, 
and of the inclination angles obtained for the precessing 
equilibrium solutions, as a function of the number of rings 
used to represent the disk (see Fig. [5}. Figure [7] shows 
how the shape of a solution near the lower stable bound- 
ary of Fig. [4] and its stability changes with the number of 
rings used to represent the disk. The transition from unsta- 
ble to stable occurs at a disk mass fraction of Md/Mth = 
0.0018,0.0013,0.0011,0.0011 for n = 15,25,45,75 rings. 
This shows that the transition has approximately converged 
when n ~ 45. 

Next, we consider the issue of scaling. We have already 
seen that the equilibrium solutions can be scaled in radius, 
mass, and precession speed according to the approximate 
(but accurate) scaling relation (|28[) : the angles 0i for all 
rings in the disk remain unchanged if the precession speed 
expressed in units of the angular frequency scales linearly 
with the disk-to-black hole mass fraction. Figure [8] shows 
disks with Md/Mbh = 0.05 consisting of equal mass rings 
with constant ratios of all ring radii relative to each other. 
Keeping the ratio 4>\fr^ ei constant, we have moved the mid- 
points of these disks to various radii. All these configurations 
are stable, and actually represent a similar warped disk (i.e. 
innermost and outermost inclinations being the same) at 
different distances from the black hole. 

Figure [5] illustrates the stability properties for these 
disks when the disk mass and precession speed are changed 
simultaneously. The boundaries of the stable solutions ob- 
tained for a fixed disk mass configuration, i.e., a configu- 
ration with constant ratios of all ring masses and radii rel- 
ative to each other, are shown by the two dashed lines in 
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Figure 8. 35-ring-disks with constant M d /M bh = 0.05 but differ- 
ent sizes. The curves from left to right show warped disk solutions 
for which the scale radius r re f is shifted to ever higher values while 

correspondingly rfi is decreased to keep rf>^J rj? ef fixed. The inner- 
most and outermost ring inclinations arc identical for all disks, 
as expected from the scaling relation, and all rescaled disks are 
found to be stable. 
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Figure 9. Schematic representation of stable warped disk solu- 
tions in the M^/Mf,^ versus <p/Q plane. The two inclined dashed 
lines enclose the stable region: the upper (lower) dashed line cor- 
responds to the maximum (minimum) stable precession frequency 
at fixed disk-to-black hole mass ratio. This diagram is for disks 
with the same surface density profile and inner and outer bound- 
aries r in , r ou t- 



the figure. All stable solutions lie between the two lines. 
These dashed lines are nearly straight, indicating that the 
minimum and maximum <fi/Q Ic f solutions at different mass 
fractions can be essentially scaled to eachother - if the same 
(scaled) equilibrium solution was the stability boundary for 
all mass ratios, the lines would reflect the scaling relation of 
eq. (|27[l respectively eq. (|28p precisely. 

Figure [9] shows that for each ratio of disk mass to black 
hole mass, there is a minimum and maximum stable preces- 
sion speed 4>/Cl, and vice-versa, (see Fig. [4] and Section f3.ll 
below). The minimum (maximum) stable mass for given <j>/Q 
is stable for all precession speeds lower (higher) than the 
original <j>/Q, until the other boundary curve is reached. For 
other disk mass configurations, the range of stable mass ra- 
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Figure 10. Inclination of a disk of 35 rings at different radii 
for Md/Mbh = 0.0016. On each curve the precession frequency <fi 
has a different value, given on the figure in terms of the orbital 
frequency at the position of the middle ring. The upper and lower 
curves show the boundaries of stable solutions. See Section 12.71 
and Fig. |9] 
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Figure 11. Variation with disk-to-black hole mass ratio of the in- 
clinations of the outermost (9 > 0) and innermost (9 < 0) rings of 
a 35-ring warped disk with r ou t/ri n = 1.44, for two different val- 
ues of the precession frequency. The figure shows how the warping 
increases with increasing disk mass fraction and with decreasing 
precession speed, and also illustrates the scaling relation of equa- 
tion ]28p . Up to numerical errors, the minimum and maximum 
values of 9 arc the same on all two curves. 



tios and the corresponding boundary lines change, but the 
qualitative behaviour remains the same. 



3 STEADILY PROCESSING, NON-LINE ARLY 
WARPED KEPLERIAN DISKS: RESULTS 

3.1 Warp Shapes and Warp Angles of Stable 
Precessing Disks 

We have already seen in the discussion of three-ring and 
n-ring systems in Section [5] that self-gravitating precessing 
disks in a Keplerian potential can be strongly warped. In 
fact, some of the disks shown in Fig.|4]are so strongly warped 
that they would obscure the central black hole from most 
lines-of-sight. 

In this section we discuss these results in more detail. 
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Figure [10] shows the warped stable equilibrium solutions ob- 
tained for a sequence of disks with varying precession fre- 
quency. These disks have constant surface density between 
fixed inner and outer radii, and a total disk-to-black hole 
mass ratio of Md = O.OOlBMt,^. The solutions shown in 
Fig. [10] are all linearly stable, according to the analysis de- 
scribed in Section [2.7l Outside the range of models bounded 
by the upper and lower curves one can find further equilib- 
ria, but these are unstable. 

By construction, these disks have a fixed line-of-nodes 
at all radii, and their shapes are given in terms of the in- 
clination angle 9 relative to the plane defined by the total 
angular momentum vector. In all cases there is a middle 
section of the disk which lies approximately in this plane, 
whereas the inner and outer parts warp away from this plane 
in opposite directions. For the most strongly warped stable 
solution in Fig. [10] the inner warp is by ~ 180 deg and the 
outer warp by ~ 120 deg. This is obtained for the lowest sta- 
ble pattern speed, in accordance with the balance between 
gravitational and Coriolis torques [see Fig. [TJ and equation 
(|19p j: the torques are weakest for the large inclinations. This 
can be seen already in the two-ring problem [see equations 
(|20[) and <|22[) ]. The least strongly warped disk solution in 
this example has inner and outer warps ~25-30deg. 

The variation of the maximum inner and outer warp 
angles with disk mass fraction is shown in Figure [TTJ for 
fixed precession frequency and radial extent of the disk. The 
curves with 8 > represent the outermost ring inclinations, 
and those with 6 < show the innermost ring inclinations 
for different precession speeds. As we have already seen in 
Figs. [J] 1101 these inclinations increase with increasing disk 
mass fraction and with decreasing precession speed. Fig. [TTJ 
also illustrates the scaling relation of equation 1)28 P- Up to 
numerical errors, the different curves can be scaled on top 
of each other. 

Some three-dimensional illustrations of warped disks 
from this family are shown in Figure 1121 From top to bot- 
tom, these plots shows warped disks with increasing ampli- 
tude of the warp, such that the disk in the bottom panel of 
Fig. [12] completely encloses the central black hole. 



3.2 Comparison with Linear Theory Solutions 

Previous work on warped and twisted disks around black 
holes has often made use of the linear approximation, in 
which the inclination angles of all parts of the warped disk 
are assumed to be small. It is therefore interesting to briefly 
consider the linear limit of our analysis above. 

For small warping angles aiij, the self-gravity torques 
simplify considerably. Because the leading term in equa- 
tion (I13al) is already 0(aij), only the 0(1) part of the Iij 
terms in this equation need to be included, while the next 
order, 0(a 2 j), can be neglected. Thus in computing Iij the 
k 2 term for ring j in equation (|10p can be approximated as: 
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and the Iij term can be integrated to give 



(30) 
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Figure 12. Three-dimensional views of several 15-ring disks sim- 
ilar to those shown in Fig. [4] with mass ratios given on the plots. 
The degree of warping increases with the disk mass fraction until 
(for disk mass fraction greater than 2 percent) the central black 
hole is completely hidden behind the warped disk. 
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Figure 13. Warp shapes for a disk with M d /M bh = 0.005 as 
obtained in the linear (dashed line), and nonlinear (solid line) 
torques regimes, for the same precession speed. The linear mode 
is scaled to the maximal amplitude for which the linear approxi- 
mation to the gravitational torques can be used. The non-linear 
disk shape is determined for the precession frequency given by the 
eigenvalue of the linear mode; it is unstable. Stable non-linear 
warps for this mass configuration have lower precession speeds 
and are more strongly warped than the solution shown. 



The mutual torque on ring i from ring j becomes, to first 
order in an, 



dVij dajj 
dctij d9i 



2GrnimjrirjIijOLij dcti 



do, 



(32) 



For a precessing equilibrium when (f>i = (f>j, 
the equation pg i = becomes to 0{9): 



Pe, 



IGrmrajTirjlij^ 

3=1 



(rf + r 2 ) 3 / 2 



(33) 

Equation (|33|l is a quadratic eigenvalue problem for the pre- 
cession frequencies. When linearized, it transforms into a 
generalized eigenvalue problem of dimensions 2n x 2n. We 
use the NAG routine F02BJF to find the eigenvalues and 
eigenvectors of equation (f33)l . 

The 2n eigenvalues constitute two distinct families in 
the frequency spectrum of the disk; the fast prograde, and 
the slow retrograde frequencies. When sorted in decreasing 
order, the first retrograde frequency has a value of zero, 
and the associated eigenvector represents a tilt of the whole 
disk by the same angle, i.e. 9i — constant. The next eigen- 
value corresponds to the warps we have discussed so far 
where the disk has one radial node (modified tilt mode, 
Hunter fc Toomrj Il969l ; ISparkei [l98i ISparke fc Casertanol 
19881 ). In the following, we restrict our discussion to linear 
warp shapes of this kind. 

In Figure [13] we show the modified tilt mode in linear 
theory of a disk with M d /M bh = 0.005, r- m = 5 and r out = 
7.2 for 40 rings. This is obtained by solving equation (|33[) 
and is shown with the dashed line. We note that in linear 
theory the warp shape can be arbitrarily scaled as long as the 
local gradient of the tilt satisfies d9/dr < 1.2/r (see Section 
(|2.2[) ; if this condition is violated, the linear approximation 
to the self-gravity torques breaks down. Therefore in Fig. [T3] 
the linear mode is scaled to its maximum possible amplitude 
such that the condition is everywhere satisfied. 
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Figure 14. Inclination of a 35-ring-disk at different radii as in 
Figure |4l but for an exponential surface density profile with scale 
length r d = 2.5. The precession speed is identical to that in Fig. 3] 
The upper and lower curves correspond to the limiting mass ratios 
for disk stability; see Section 12.71 



As mentioned above, the precession frequency of this 
mode is the first nontrivial eigenvalue in the retrograde fam- 
ily, and here it has a value of (/>/fi ro f = —0.0103 when nor- 
malized to the rotation frequency of the reference ring. For 
this frequency, we then solve equation (|5a|l to obtain the 
nonlinear warp shape shown by the solid line in Figure [T3l 
The larger curvature of the non-linear warp near the inner 
and outer boundaries of the disk, with respect to the scaled 
linear mode, shows that the linear approximation overesti- 
mates the torques in these parts of the disk. 

However, the main difference between linear modes and 
non-linear warps is that, for a given mass distribution of the 
disk (surface density profile, inner and outer boundaries), 
the precession frequency and shape of the modified tilt mode 
in linear theory is uniquely determined, whereas non-linear 
equilibrium warp solutions may exist for a range of preces- 
sion frequencies and warp shapes or, e.g., for extended disks, 
may not exist at all. For the case shown in Fig. 1131 non- linear 
warped equilibria are found for precession frequencies in the 
range 0/f2 rcf = -8.87 x 10~ 5 — ► -1.16 x 10" 2 and are stable 
in the range 0/n rcf = -9.24 x 10 -5 -> -9.63 x 10" 3 . The 
particular non-linear warp shape obtained for the frequency 
of the linear mode and shown in Fig. [13] is unstable. 

Alternatively, the warp shape may be parametrized by 
the inclination of the outermost ring, say, 9 n . Linear the- 
ory warps can in the previous example be considered valid 
up to 9 n ~ 10° , and have all the same precession frequency. 
Non-linear warp modes are found in the range 6 n — 19.1° — > 
130°, and are stable in the range 6„ = 20.9° -> 129°. They 
are disjunct in 9 n from the linear modes, and their preces- 
sion speed decreases with 9 n according to the balance of 
gravitational and Coriolis torques. 



3.3 Dependence on Surface Density Profile and 
Radial Extent 

The warped disks presented in Figs. [4] [TT] have constant 
surface density. For comparison, Figure [T?] shows the warp- 
ing of an exponential disk, with surface density E(r) = 
En ex P(— r/rd) where En denotes the central density and r d 
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Figure 15. Stable warped equilibria for 35-ring-disks with con- 
stant surface density but varying radial extent, and with Md = 
0.01M bh . The solid and dashed lines demarcate the boundaries of 
the region of stable solutions. 
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Figure 16. Maximum radial extent of stable warped disk equi- 
libria r ut /fin, as a function of disk-to- black hole mass ratio. The 
dashed line represents the result for a disk of 15 rings with con- 
stant surface density, the solid line is for an exponential disk with 
scale-length = 5.2, and the dot-dashed line is for an exponen- 
tial disk with scale-length = 2.5. The inset shows results ob- 
tained for the constant surface density disks when instead approx- 
imated by 45 rings, near the mass ratio with maximum r out /ri n . 
Typical ratios are r ou t/r; n =2-4. 



is the scale length, chosen to be 2.5 units in this example. 
The other parameters (relative ring radii, precession speed) 
are identical to those used in Fig. [4] The basic warp shapes 
are similar as for constant surface density, but the maximum 
outer warp angles are slightly larger. The range of stable disk 
masses is also comparable to that for the constant surface 
density disk (for the same precession speed); see the curves 
showing the boundaries of stability in Figs. [4] 1141 

Because the condition for a warped equilibrium is that 
the Coriolis and gravitational torques balance, clearly not 
only the mass fraction and mass distribution, but also the 
radial extent of the disk must be important for determin- 
ing the warp shape and its stability. To investigate this we 
compiled a set of precessing equilibria with varying radius 
scaling factor re, as follows (see also Section r2.5|l . After fixing 
the radius of the middle ring of the disk, r m jd, we determine 



J = i 
3 = i ~ if 



n + l 
2 



n odd 
ii even 



1 = 1,2, 



,n (34) 



where n is the number of rings. For illustration, we consider 
a family of disk models with the same disk-to-black hole 
mass fraction Md = 0.01M(,h, each with its own constant 
surface density given by Md and k. All disks are made of 
n — 35 rings, and the middle ring radius is set to r m id = 6 
units. 

Figure [Tol shows precessing equilibria for such disks for 
different ft. The upper and lower curves show the two disk 
shapes that bound the stable range of solutions in terms 
of the /t-factor. In the case where the rings have minimum 
possible separation from each other, the inner ring has a 
radius of n = 5.7 units, and the outer ring has r n = 6.3 
units. On the other hand, for the most extended stable disk 
in this family, n = 3.9, and r n — 9.1. When the extent 
of the disk is increased, a slight decrease in the warping is 
observed in Fig. [15] This is due to the fact that the torque 
from a ring of constant mass decreases with distance to the 
ring, cf. equation (|13a|l . 

Figure [16] shows the radial extent of the disk r out /n n 
for which stable warped equilibria can be found, for dif- 
ferent surface density profiles and as a function of disk-to- 
black hole mass ratio. The most important result of these 
calculations is that stable non-linear warps can be main- 
tained only for disks with inner and outer boundaries, for 
which r out /r in — 2 -4. T his is reminiscent of the result of 
iHunter fc Toomrei | | 19691 ) that in linear theory only trun- 
cated disks permit long-lasting bending modes. 



3.4 Time-Evolution of Ring Systems 

In this section we consider the explicit time-evolution of a 
precessing system of self-gravitating rings in a massive black 
hole potential. By integrating the equations of motion, equa- 
tions (J3j)- (J6jl , starting from initial conditions corresponding 
to one of the precessing disk solutions found earlier, we can 
check the stability of this solution directly and compare with 
the linear stability analysis. 

In these integrations, we use disks of 20 equal mass 
rings, equally spaced in radius. The ratio of the outermost 
ring radius to that of the innermost ring is 1.44. The initial 8i 
are obtained from precessing equilibrium solutions; all rings 
have the same line-of-nodes, i.e., the same initial <f>i. The 
equilibrium precession speed is given by <fi/tt YO f — —0.0021. 

In the following figures, symbols starting from the outer 
circle show the variation of inclination 9 with ring radius, 
where the ring radii are shown as distances from the center 
of the plot, with scale shown on the lower right. The symbols 
starting from the inner circle show how the azimuthal angle 
4> changes with the ring radius; for this part of the plot, 
the ring radii are scaled down to the half of their values to 
make the figure more easily readable. The elapsed time of 
the integration is shown on top of the figures, in terms of the 
number of orbital periods at the position of the outermost 
ring n where </>/ft n = -0.0027. 

Figure [T7] shows the time evolution of a disk of 20 rings 
with Md = 0.05Af(,h. The disk stays in equilibrium for 12 
orbital periods, consistent with its linear stability. Figure [T8l 



12 A. Ulubay-Siddiki, 0. Gerhard & M. Arnaboldi 



starting configuration 4 orbital periods starting configuration 8 orbital periods 

• e • e e 



1 2 3 4 5 6 7 



1 2 3 4 5 6 7 



1 2 3 4 5 6 7 



1 2 3 4 5 6 7 



8 orbital periods 



12 orbital periods 



1 orbital periods 



1 6 orbital periods 



1 2 3 4 5 6 7 



1 2 3 4 5 6 7 



1 2 3 4 5 6 7 



1 2 3 4 5 6 7 



Figure 17. Time-evolution of a disk of 20 rings with total mass 
M d = 0.05M bh . The disk is stable. 



Figure 18. Time evolution of a disk of 20 rings with total mass 
M d = O.lMi,^. The disk precesses as a unit for nearly 8 orbital 
periods, but then it starts to break up, so is unstable. 



shows the evolution of disk of 20 rings with M d = O.lMth- 
This disk precesses as a unit for 8 orbital periods, but then 
it starts to break into parts, hence the disk is unstable, as 
also predicted by linear stability analysis. 

To strengthen the agreement between linear stability 
and time evolution results, we integrate two of these ring sys- 
tems for longer. Figures [19] and [20] show the time evolution of 
two disks with masses M d = 0.005M b/i and M d = 0.02Mi, h , 
with 15 logarithmically spaced rings according to equa- 
tion (|34[) . over 50 orbital periods. In both cases the disks 
are stable, as expected from the linear stability analysis. 



starting configuration 



50 orbital periods 



. ♦ 



4 DISCUSSION 

4.1 Theoretical Issues 

In this work we have considered warped disks around black 
holes for which the only acting force is gravity and the 
disk is approximated as a nested sequence of circular rings. 
We have focussed on non-linearly warped, steadily process- 
ing disk configurations, contrary to most previous work in 
which small amplitude warps were considered, often of a 
transient nature. We have found that stable, steadily pro- 
cessing, highly warped disks can be constructed, albeit only 
over a limited radial range, such that the typical ratio of the 
outer to the inner boundary radius is ~ 2-4. 

In one illustrative case, we have compared with a lin- 
ear theory warped disk. For a given disk mass configuration, 
the precession frequency of the linear, modified tilt mode is 
given as an eigenvalue, and the shape can be scaled up to 
the amplitude where the validity of the linear approximation 
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Figure 19. Time-evolution of a disk of 15 rings with mass frac- 
tion Md/M^h = 0.005. The disk is followed for 50 orbital periods 
and is stable. 

to the gravitational torques breaks down. The corresponding 
non-linear warp with the same precession frequency is unsta- 
ble. Stable non-linear warps for the same mass configuration 
exist for a disjunct range of precession speeds which are all 
slower than that of the linear mode. Their warp angles in- 
crease with decreasing precession speed, and the non-linear 
solutions are more strongly warped than the linear mode at 
the maximum scaling. 

These warped disks obey a scaling relation in the sense 
that (i) they can be scaled to an arbitrary radius r, provided 
the precession speed is scaled to the circular frequency fi(r), 
and (ii) they can be scaled in mass, provided the ratio of 
precession frequency to Q(r) keeps in line with the ratio of 
disk mass to black hole mass. 
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' e 



next step is to add the gravitational torques to these equa- 
tions and study the evolution of viscous, self-gravitating, 
non-linearly warped disks; this is work in progress. 



01234567 01234567 
Figure 20. Time-evolution of a disk of 15 rings with mass frac- 
tion of M d /M bh = 0.02. The evolution is followed for 50 orbital 
periods and the disk is stable. 
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Figure 21. Inclination of a disk of constant surface density at 
different radii, under the self-gravity torques and the quadrupolc 
torques of a surrounding star cluster. This is parametrized by the 
ratio TjMt/Mj, where r\ mea sures the flattening of the potential 
(equation (2) of[Sparkc 198(|) and M t /M<i is the ratio of the mass 
of the star cluster to the mass of the disk, inside the outermost 
ring radius. 



In constructing these solutions, we have neglected the 
background potential generated by the surrounding nuclear 
star cluster, whose quadruple moment will often be impor- 
tant on scales of ~ 0.1 pc. Figure [5T] shows steadily pro- 
cessing warped disk inclinations for one case including the 
background potential. For the parameters chosen, the solu- 
tions are qualitatively similar to those discussed earlier. 

Stability was tested with respect to perturbations of the 
ring parameters, that is, the orbits of gas and stars were as- 
sumed to remain circular. We did not investigate instabilities 
by which the disk would become eccentric or lop-sided. An- 
swering the question whether such instabilities are relevant 
for the warped disks considered here requires differe nt tech- 
nique s and must remain for future work (see, e.g., iToumal 
12002ft . 

Neglecting gas pressure and viscosity for our warped 
disk solutions is justified if the disks are cold and the viscous 
time-scale is much longer than the precession time-scale. 
IPringld (QUI) has devised a system of equations for the 
evolution of the surface density and local angular momen- 
tum vector of a non-linearly warped, viscous disk. A logical 



4.2 Origin of Warped Disks 

An important question is whether, and if so how, the warped 
nuclear disks we have considered can be set up in nature. 
Infall of gas clouds on inclined orbits has been discussed in 
the context of observations of the Galactic center (see next 
section) as a possi ble model for generating a warped disk in 
the central parsec (|Hobbs fe Navakshin|[2TM ). If the poten- 
tial of the nuclear star cluster is important, accretion of gas 
onto a plane inclined relative to its principal plane may lead 
to a warped disk. The combined quadrupole moment of the 
gas disk itself and of the background cluster potential would 
cause the orbits to precess and the disk to become warped. In 
both cases, the accreting gaseous material with misaligned 
angular momenta will not directly end up in a warped disk 
with the right density structure for steady state precession. 
However, the disk may settle into a warp mode if the energy 
associated with the transie nt response can be transported 
outw ards by bending waves (|Toomrel ll983l ; lHofner fe Sparkel 
1 1994T ) . or in the c ase of gaseous disks, if it can be dissipated 
(see discussion in lPapaloizou et al.lll99Sf ); this remains to be 

inve stigated. 

ICaproni et all (|2006t ) discuss four warping mechanisms 
for extragalactic accretion disks: tidal, irradiative, magnetic, 
and Bardeen-Petterson. If a planar disk has beco me warped 
by the radiation pres sure instability discussed by iPetterson 
lll97^ ; |pringie]||l996t ) or through magnetic instabilities (jLail 
120031 ) . the gravitational torques might start to dominate 
once the source of the initial warping disa ppears. Highly 
warped disks have been reported before by IPringld (|l997T ) 
in the context of the radiation pressure instability. We have 
done some simple time evolution calculations to show that 
initially highly warped disks often do not dissolve through 
self- gravity precession; the torques then cause wobbling but 
not break-up of the disk. The role of self-gravity in such 
models would be to ensure the long term persistence of the 
warp. Future work along the lines discussed at the end of 
the last subsection may be able to clarify whether this is 
feasible. 



4.3 Warped Disks in Galactic Nuclei 

Warped disks around central black holes have been inferred 
through observations of water maser emission in several 
nearby active galaxies such as NGC 4258, NGC 1068, and 
the Circinus galaxy. The maser dis ks in these galaxies ex - 
tend radially b etween 0.16 — 0.28 pc fllerrnstein et al.lll999l ) , 
0.65 - 0.11 pc dGreenhill fe Gwinn|[l997l ) . and 0.11 - 0.4 pc 
(|Greenhill et al.ll2003bh T respectively. The most widely stud- 
ied of these maser disks is in NGC 4258, where from the 
near-Keplerian rotation curve of the high-velocity masers 
the black hole mass is deduced to be 3.8 x 10 7 Mo, and 
the dyn amical upper limit to the mass of the disk is < 
1O 6 M (|Herrnstein et al.ll2005t) . Stationary, power-law ac- 
cretion disk models constrained by theory and observations 
have mass fractions 10 -4 — 10~ 3 of the central black hole, 
in which case the gravitational and viscous torques are 
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comparable (|Caproni et alj|2007l : iMartml [20(51 ) . Several ex- 
planations have been suggest ed for the observed warp in 
the disk (jCaproni et alJ l2006t >. In one model the warp is 
caused by a binary co mpanion orbiting outside the disk 
l|Papaloizou et al1ll998l ) ; this would need a mass comparable 
to that of the disk. A se cond possibility is radiation pressur e 
fro m the central source ( Pring lc 1996; Maloi iev et al.|[l996h . 
but ICaproni et al.l <|200C ) analysing several AGN disks find 
that these are stable against radiation warping. The most 
favour ed explanation for th e warp is the B ardeen-Petterson 
effect (|Caproni et alJ 120071 ; iMartinl 120081 ) but to reach a 
steady state the disk must be very long-lived. Gravitational 
torques have so far been mostly neglected; our results sug- 
gest that it may be worthwhile to consider models including 
both the gravity from the disk and possibly the quadrupole 
moment of the stellar cusp. 

In the Galactic center, near-infrared observations have 
identified one or possibly two disks of young stars at a dis- 
tance of ~ 0.04 to . 4 pc from the central b lack hole SgrA* 



dGenzel et al 



iBartko et al.l 



120031 ; iPaumard et~aH 120061 ; iLu et all 120091 ; 
2009t ). These stellar disks are highly inclined 
both with respect to the Galactic plane, and with respect 
to each other. The total mass in the disks, as inferred from 
stella r number counts, is around 10 4 Mq (|Paumard et al.l 
120061 ). This is a non-n egligible fraction of the mass of SgrA* , 
M bh ~ 4 x 10 6 M Q dGenzel et al.ll2000l; IChez et al.ll2005l) . 
The recent analysis of IBartko et al.l (|2009l ) shows that the 
clockwise rotating disk is warped, with angular momentum 
direction slewing over ~ 60° from the inner to the outer 
stars. We consider the precession of the warped disk in the 
Galactic center elsewhere. 

Warped disks could al so have impor tant implications 
for the unification of AGN (|Phinnevlll989l ). The unification 
theories rely on the obscuration along some lines-of-sight of 
the radiation from the central source by intervening mat- 
ter. While this obscuring matter is usually depicted as a 
doughnut-like torus, an alternative possibility is that it could 
have the shape of a flared or warped disk. The highly warped 
solutions discussed above in principle provide the geome- 
try to obscure the central engine from most lines-of-sight. 
The obscuring medi um required for these unification scenar- 
ios must be clumpy ([Nenkova et al.l 20021), perh a ps suggest- 
ing fr agmentation of the disk ( Goodman! 120031 ) . iNavakshinl 
(|2005l ) studied the evolution of a highly inclined warped disk, 
where he showed that the disk indeed can conceal the cen- 
tral object for most of its lifetime. In the nonlinear regime, 
warped disks can obscure a significant part of the solid angle 
of the source (s ee Figure [T2l in Section [3] above). Recently, 
IWu et aTl (|2008h showed that because the outer parts of a 
warped disk receive a larger fraction of the central emission, 
the line ratios of the reprocessed Balmer emission lines can 
be successfully predicted by a warped disk model. 



5 SUMMARY AND CONCLUSIONS 

In this paper we have investigated non-linearly warped disk 
solutions around black holes for which the only acting force 
is gravity. We used a simple model in which the disk is ap- 
proximated as a nested sequence of circular rings. We have 
shown that with these approximations stable, steadily pro- 
cessing, highly warped disks can be constructed. 



These disks have a common line-of-nodes for all rings. 
In all cases there is a middle section of the disk which lies 
approximately in this plane, whereas the inner and outer 
parts warp away from this plane in opposite directions. The 
warp angles of these solutions can be very large, up to ~ 
±120 deg, but they extend only over a limited radial range, 
such that the typical ratio of the outer to the inner boundary 
radius is ~ 2-4. Such precessing equilibria exist for a wide 
range of disk-to- black hole mass ratios Md/Mbh, including 
quite massive disks. 

The stability of these precessing disks was determined 
using linear perturbation theory and, in a few cases, con- 
firmed by numerical integration of the equations of motion. 
We found that over most of the parameter range investi- 
gated, the precessing equilibria are stable, but some are un- 
stable. 

These disks obey a scaling relation: they can be scaled 
to arbitrary radii r, provided the precession speed is scaled 
to the circular frequency fl(r), and they can be scaled in 
mass, provided the ratio of precession frequency to Q(r) is 
changed, in good approximation, proportionally to the ratio 
of disk mass to black hole mass. 

The main result of this study is that persistent forc- 
ing of the disk other than by its own self-gravity is not 
necessarily required for maintaining a non-linearly warped 
disk in a Keplerian potential. Further work combining self- 
gravity with gas physics etc. will show whether these self- 
gravitating warped disk solutions help to understand the 
observed warped disks in galactic nuclei. 
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